Breaking conjugate pairing in thermostatted billiards by magnetic field 
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We demonstrate that in the thermostatted three-dimensional Lorentz gas the symmetry of the Lya- 
punov spectrum can be broken by adding to the system an external magnetic field not perpendicular 
to the electric field. For perpendicular field vectors, there is a Hamiltonian reformulation of the dy- 
namics and the conjugate pairing rule still holds. This indicates that symmetric Lyapunov spectra 
has nothing to do with time reversal symmetry or reversibility; instead, it seems to be related to 
the existence of a Hamiltonian connection. 
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I. INTRODUCTION 

Thermostatted dynamical systems have raised consid- 
erable interest recently as a testing ground for ideas in 
nonequilibrium statistical mechanics |EJ. In particular, 
questions concerning the role played by chaotic dynam- 
ics in the appearance of nonequilibrium stationary states 
in dissipative systems have been in the focus of research 
activities [0. One of the most remarkable features of 
these models is that they are dissipative and time re- 
versal symmetric at the same time. Some, but not all, 
thermostatted systems have another interesting common 
property known as the conjugate pairing rule (CPR): the 
Lyapunov exponents of the system form pairs summing 
up to the same (negative) value || . It has a practical rel- 
evance too: with CPR, one pair of Lyapunov exponents 
can be used to determine the sum of all the exponents, 
which is known to be connected to the transport prop- 
erties of the system j| . CPR is trivially present in con- 
servative Hamiltonian systems (the sum being zero due 
to simplecticity); however, there is no obvious reason to 
expect anything similar in dissipative systems. In fact, 
CPR in thermostatted systems was first discovered by 
numerical studies [|| . 

The simplest system in which CPR can be checked 
is the three-dimensional (periodic) Lorentz gas (3DLG): 
due to its three degrees of freedom, it has four nontrivial 
Lyapunov exponents. Dcttmann et at have shown nu- 
merically |H that the 3DLG with external electric field 
and Gaussian isokinetic (GIK) thermostat exhibits con- 
jugate pairing; later, this has been proven analytically 
for conservative forces and hard- wall scatterers H. It 
has also been demonstrated [JtJ that this system can be 
connected to a Hamiltonian dynamics. 

In this paper, we check the effect of an external mag- 
netic field on the validity of CPR in the GIK thermostat- 
ted cubic lattice 3DLG. In particular, we will focus on 
two features possibly related to CPR: reversibility (an 
extension of time reversal symmetry) and the existence 
of a Hamiltonian formulation. Both can be controlled by 
the direction of the magnetic field with respect to that of 
the electric field and the lattice. Our numerical results 
show that CPR is not affected by breaking reversibility, 



and it also holds for cases with perpendicular electric and 
magnetic field vectors for which there is a connection to 
Hamiltonian dynamics. However, CPR breaks down for 
nonperpendicular fields, i.e. in the case when no Hamil- 
tonian connection has been found. 

In Sec. |J, the equations of motion for billiards with 
GIK thermostat in magnetic field are presented, together 
with a discussion of reversibility and the Hamiltonian 
connection for perpendicular fields. The numerical re- 
sults for the 3DLG and our conclusions are presented in 
Sec. Ill and [V, respectively. 



II. 



THERMOSTATTED BILLIARDS IN 
MAGNETIC FIELD 

A. The dynamics 



The kinetic energy of a particle moving under the in- 
fluence of external fields can be kept constant by adding 
a special frictionlike force to the system. Since this force 
can be deduced from Gauss's principle of least constraint, 
and it is kinetic energy that is kept constant, the tech- 
nique is called Gaussian isokinetic (GIK) thermostat . 
In billiards, this is equivalent of particle momentum p 
changing only in direction, but not in magnitude p, dur- 
ing the "free" flights between collisions with the hard-wall 
boundaries. Choosing the unit of mass to be the mass of 
the particle, the corresponding equations of motion are 



P, P 



ap 



(1) 



where q = (x, y. z) is the position of the particle, F e 
stands for the external forces, while the GIK thermostat 
corresponds to the choice 



a 



P 



(2) 



For simplicity, we will choose length and time units in our 
studies so that p = 1, but care must be taken when sub- 
stituting 1 for p 2 in terms like a above, especially in the 
derivation of tangent space equations for the calculation 
of Lyapunov exponents. 
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In our model, the external force F e contains the (con- 
stant) electric and magnetic fields E and B: 

F e = E + p x B (3) 

(we have defined the unit of electric charge to be that of 
the particle) . The full dynamics also includes the secular 
collisions with the hard-wall boundaries, changing the 
momentum pi to p/ instantaneously: 

p / = (/-2non)p ! (4) 

where / is the (3 x 3) identity matrix, n is the normal vec- 
tor of the boundary at the collision point and 'o' denotes 
the diadic product. In our 3DLG, the scatterers are hard 
spheres of radius R, arranged into a regular cubic lattice 
with distance d between the centers of nearest neighbor 
scatterers. For simplicity, we choose the length scale so 
that R = 1. 



B. Reversibility 

Without magnetic field, Eqs. (]l|) and (||) ensure time 
reversal symmetry for the dynamics, which means that 
for each solution T + (t) = (q(£),p(i)) T there exists an- 
other one tracing the same path backward in time: 

L(t) = (qH),-pH) T (5) 

The pairing of solutions by time reversal symmetry is im- 
portant in these models: it is used e.g. in showing that 
the average current flows in the direction of the external 
electric field This symmetry cannot hold if B ^ 0, 
but the more general property of reversibility Q may still 
be true, depending on the particular choice of E and B. 
Reversibility means the existence of a transformation G 
in phase space which is an involution (i.e., G 2 is the iden- 
tity) mapping each solution r+(t) to another one T_(i) 
in the following manner: 

T-(t)=GT+{-t). (6) 

In terms of the the phase space flow defined by T(t) = 
<fi T(0), this requirement can be written as 

GtfG = jr* , (7) 

i.e., bracketing the flow by G "reverses the direction of 
time" . 

Ordinary time reversal symmetry is equivalent to 
G = Gq just flippping the direction of the momentum: 
Co(q, p) = (q, — p) T - For B ^ 0, the flow can be re- 
versed by the transformation Gb — MGq where M is a 
mirroring of q and p with respect to the plane contain- 
ing E and B (the proof of this statement is left to the 
Appendix). In the Lorentz gas, reversibility of the full 
dynamics also requires that the invariant plane of M be 
a symmetry plane of the lattice too. This gives us an easy 
way to control reversibility in the Lorentz gas: choosing 
directions for E and B in a symmetry plane of the lat- 
tice leads to reversible dynamics, otherwise we have no 
reversibility. 



C. Hamiltonian formalism 

A nontrivial result for GIK thermostatted systems 
without magnetic field is that a Hamiltonian formulation 
of the dynamics exists provided the force F e is the gra- 
dient of a scalar field — $(q) 0. Then there is a Hamil- 
tonian H (Q, P) so that the GIK equations of motion for 
the physical variables q and p can be obtained from the 
canonical equations of motion for Q and P through a 
suitable coordinate transformation. It is straightforward 
to check that the Hamiltonian H(Q, P) = ^[e*P 2 - e - *] 
has canonical equations leading to Eq. (0) if one assumes 
the transformations q = Q and p = e*P. However, 
it is important to stress that this connection holds only 
if we make explicit use of the constraint p = 1 and its 
equivalent H = in the GIK and canonical equations, 
respectively. 

The extension of the Hamiltonian formulation to cases 
with B ^ is not as obvious as for conservative systems 
because of the factor e* in front of P 2 in the "kinetic 
energy" term of the Hamiltonian. Nevertheless, we may 
still follow a similar route by defining ^(q) through E = 
— as usual and replacing P by P — a(q) in H, where 
the vector a(q) is connected to the magnetic field. This 
leads to the Hamiltonian 

ff B (Q,P) = i[e*(P-a) 2 - e -*]. (8) 

A lengthy but straightforward calculation shows jit)) that 
the canonical equations for Hb — can be connected to 
the GIK equations of motion for p = 1 by the transfor- 
mation 

q=Q, p = e*(P-a) (9) 

if we assume the following relationship between B and a: 

B = e*rot a. (10) 

Note that this is an extension of the usual relationship 
B = rot a for the GIK thermostat. However, due to 
the presence of e* in Eq. (|l0|), we do not necessarily 
have a solution a for arbitrary E and B. Indeed, since B 
must satisfy Maxwell's equation div B = 0, this condition 
leads to the restriction EB = 0. Therefore, we can use 
the Hamiltonian formulation given above only in the case 
when B is perpendicular to E. 

III. NUMERICAL RESULTS 

We have calculated numerically the Ljapunov spec- 
trum of the GIK thermostatted 3DLG with constant 
electric and magnetic fields. The Lyapunov exponents 
Ai > A2 > • • • > A6 can be measured by following the 
evolution of a full set of linearly independent tangent 
space vectors along a very long trajectory and applying 
repeated reorthogonalization and rescaling to them; see 
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FIG. 1. Time evolution of the measured Lyapunov ex- 
ponents in the 3DLG. The scatterers are balls with ra- 
dius R = 1 arranged in a cubic lattice with lattice con- 
stant d — 2.3. The external fields are E = (0.3, 0, 0) and 
B = (0, 0.5 cos <f>, 0.5 sin <f>) with <j> — n/20. The arrows show 
the long-time values of the exponents (t — 10 s ); A3 and A4 
converge to as expected. 



Rcf. for a detailed description of this method. The 
effect of collisions were taken into account by the formula 
presented in Ref. [Q. In all cases studied we have ob- 
tained finite time exponents converging to their infinite 
time limits as in the example plotted in Fig. [j]. The fluc- 
tuations in the measured values typically tend to zero as 
1 / vN, where N is the number of collisions, so for reliable 
results we needed very long runs with N = 10 7 collisions 
or more. The data also show that the largest Ljapunov 
exponent is positive, i.e. the motion is chaotic, and that 
two of the exponents are zero as expected. 

We have choosen the coordinate axes x, y and z aligned 
with the lattice axes. Through the directions of the field 
vectors, we can have reversible or non-reversible dynam- 
ics in our model, with or without a Hamiltonian repre- 
sentation, independently. In the simulations, we fixed E 
along the x axis, so that it lies in the symmetry planes 
y = and z = 0, and controlled the above properties by 
chosing the direction of B accordingly. In particular, the 
dynamics is reversible e.g. for B y = 0; meanwhile, there 



exists a Hamiltonian formulation as given in Sec. II C for 
B x = 0. 

Figure [j] shows the results of a simulation for B = 
(0, B cos 4>, B sin </>) with <fi = 7r/20, i.e., for perpendicu- 
lar fields and without reversibility. In Fig. ^|, we plotted 
the sums of the pairs Ai + A6 and A2 + A5. They both 
converge rapidly to the same value, thus CPR seems to 
hold in this case. It is also worth noting that the differ- 
ence between the two sums disappears much faster than 
the fluctuations in the individual exponents. We have 
obtained similar results for other values of the angle (j>, 
including reversible flows (e.g. cf> — 0). These results 
demonstrate that reversibility is not needed for CPR to 
hold. 

In the second type of simulations we have choosen 




FIG. 2. Time evolution of the sums Ai + Ae (solid line) and 
A2 + A5 (dashed line), as obtained from the data in Fig. 
The two graphs are practically indistinguishable for t > 5000. 
The arrow shows the asymptotic value of the sums (t = 10 s ). 



FIG. 3. Same as Fig. 2 for nonperpendicular fields; 
B = (0.5 sin 4>, 0, 0.5 cos <f>) with <f> = 77r/20. Notice the change 
of scales with respect to Fig. 2. 



B = (B sin <fi, 0, B cos </>), so that for (f> ^ the two field 
vectors are not perp endicular and the Hamiltonian for- 
mulation of Sec. II C does not apply. The numerical Lya- 
punov spectrum looks qualitatively the same as in Fig. [l], 
but the sums of the two pairs seem to converge to differ- 
ent values as shown in Fig. || for the angle <j) = 77r/20. 
In other words, CPR is broken in this case; other val- 
ues of <p 7^ have lead to similar results. The differ- 
ence between Figs. ^| and ^| is very clear: the quantity 
A = Ai + A 6 — (A 2 + A5) converges to zero quite fast if 
CPR holds, while it stays definitely positive in the case 
without CPR. 



IV. CONCLUSIONS 

We have demonstrated that in the GIK thermostatted 
3DLG the CPR can be broken by an external magnetic 
field which is not perpendicular to the electric field. For 
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perpendicular fields, however, CPR holds, and the con- 
vergence of the pair sums to each other seems to be much 
faster than that of the individual exponents indicating 
that CPR is valid for all times in these cases just as in 
the 3DLG without magnetic field ||. This phenomenon 
is called strong CPR. The perpendicular cases are also 
characterized by the existence of a Hamiltonian formu- 
lation. There exist other nontrivial examples for sys- 
tems with strong CPR and a Hamiltonian formulation, 
too: e.g., the Gaussian isoenergetic thermostat with a 
special interparticle potential jl3| or the ideal Shod gas 
0,0. These examples suggest that there may be a di- 
rect connection between strong CPR and the existence of 
a Hamiltonian formulation. We will examine this ques- 
tion in a separate paper [(Tof . Although we are not aware 
of any counterexamples, the question concerning the ex- 
istence of systems with strong CPR but without a Hamil- 
tonian formulation is still open. 

Our results also show that time reversal symmetry, or 
reversibility in general, is not needed for CPR to hold. 
Indeed, one of the first examples for CPR in dissipative 
systems has been a Hamiltonian system with a constant 
viscous damping |Te| ] which has no time reversal symme- 
try. 
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APPENDIX 

We show that the flow defined by Eqs. is re- 

versible with res pect to the transformation Gb = MGq 
as given m Sec. |lIB|. Equation @) can be rewritten for 
the time derivative F of the flow as GFG = —F. From 
Eq. (0) one can see that F(q, p) = (p, f ) T , whith f (q, p) 
given by the expression for p. Now we can write that 



GbFGb{q, p) = G B F(Mq,-Mp) 

= G B (-Mp,f(Mq,-Mp)) 
= (-p,-Mf(Mq,-Mp)) T 
= (-p,-f(q,p) T 

The last equality gives us the condition 

f (q, p) = Mf (Mq, -Mp) 



(11) 



for the force acting on the particle. 

Since f consists of the two parts of the Lorentz force 
and the thermostat, we can check these terms separately. 
For the electric field this means that E = ME, i.e., E 
must be in the invariant plane of M. For the term p x B, 



the right hand side of Eq. ( O ) reads as M(-Mp x B) = 
-M(p|| x B - p ± x B) = -M(p|| x B) + M(p ± x B), 
where py and p± denote the components of p parallel and 
perpendicular to the invariant plane of M, respectively. 
If B is in this plane, then — M(pii x B) = p|| x B and 
M(pj_ x B) = p^ x B, so the magnetic part of the Lorentz 
force also satisfies Eq. (|TTJ) . As for the thermostatting 
force, M((E • Mp)Mp) = ((E • Mp)M 2 p = (E • p)p 
also holds if E = ME. Thus the flow is reversed by 
Gb = MGo if the field vectors E and B are invariant 
under M. 
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